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Page  7 '  The 

incorrectly, 
as  follows: 


important  derivatives  and  are  obtained 

The  corrected  values  for  these  derivatives  should  read 


JtA  e^Cbc)  -  m2h(bc) 
2  m2  —  1 


(4a) 


,  _  itA  EM  (BC)  -  M2h(BC) 

16  M2  -  1 


(5&) 


The  same  correction  should  be  made  for  these  derivatives  in  table  I 
in  the  column  "Principal  body  axes."  In  the  column  “Stability  axes" 
of  table  I  the  corrected  expressions  are 


nA  E"(BC)  -  MgH(BC) 
2  M2  _  1 


.4(1  +  8  ^SV'(BC)  -tfH(BC) 
^  16  \  0  /  M2  -  1 


The  error  in  the  derivation  consists  in  the  assumption,  carried 
over  from  reference  1,  that  the  surface  potential,  equation  (3),  1b 
not  altered  by  a  small  normal  acceleration.  This  assumption  is  true 
for  the  narrow  triangles  treated  in  reference  1  but  fails  for  the 
general  triangles  treated  in  the  present  paper.  The  assumption  is 

equivalent  to  the  neglect  of  the  time  dependency  terms  ^  0^ 

a2 

and  Mt  (where  t  represents  time  and  a  is  the  speed  of  sound) 
8, 

in  the  linearized  partial  differential  equation  for  unsteady  motions: 
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b20xx  -  0yy  ~  0zz  +  =1  ^  +  i  kt  =  0  (a) 

El  8. 

which  leaves  only  the  steady— state  (Prandtl-Glauert)  equation 


(b) 


The  correct  potential  to  replace  equation  (3)  must  satisfy  both 
equation  (a)  and  the  boundary  condition  on  the  wing 


=  —  aVt 


(c) 


In  an  unpublished  paper,  Mr.  Clifford  S.  Gardner  has,  in  effect,  shown 
that  a  suitable  solution  is 


i 

m 

a 


(d) 


where  is  the  steady— state  potential  corresponding  to  a  unit 
pitching  velocity  about  the  y-axis  and  x  is  the  steady— state 
potential  corresponding  to  unit  angle  of  attack.  Both  \|/  and  X 
satisfy  equations  of  the  form  (b).  That  equation  (d)  is  a  solution 
can  be  verified  by  direct  substitution  into  equations  (a)  and  (c). 
Thus,  Gardner  has  shown  that  the  tjme— dependent  potential  for  an  angle 
of  attack  dt  may  be  compounded  of  two  time— free,  or  steady-state, 
potentials,  one  for  a  constant  angle  of  attack  and  the  other  for 
steady  pitching. 


The  lift  distribution  at  time  t  =  0  for  the  angle  of  attack  at 
is  obtained  from  the  potential  by 


AP  =  2p(V0x  +  0t) 


=  2pV<x(M|  |  fex.-  X 


B' 


Vb2  x 

m2(ap) 


VBC 


VB£ 


_  M_£(ap)  _  AP*! 

q=l  j  v  'a=l  ^  | 


(e) 
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where 

(AP)g_i  lift  distribution  for  unit  pitching  velocity 

(AP)  1  lift  distribution  for  unit  angle  of  attack 

AP*  value  of  AP  per  unit  a  used  in  deriving  the 

incorrect  equations  (4)  and  (5) 

Integration  of  equation  (c)  over  the  plan  form  to  obtain  the  lift  and 
moment  and  reduction  to  coefficient  form  yields 


+  2 


(f) 


.2/  AP 


plan 

form 


ipv2 


dx  dy 


where  the  *  designates  the  incorrect  values  in  equations  (4)  and  (5) 
respectively,  and  the  y-axis  is  taken  through  the  center  of  gravity. 
Values  of  Ct  ,  Cm  ,  and  Cm  ,  may  be  obtained  from  table  I 
*  1  ^  *  u 

and  and  *  from  equations  (4)  and  (5).  The 


quantity 


(  \  is  obtained  by  setting  a  =  1  and  a  = 

(kL  b 


c  +  X 


in  equation  (l) .  Substitution  and  integration  then  yields  the 
corrected  values  for  Cj^  and  (equations  (4a)  and  (5a))  as  set 

forth  at  the  beginning  of  this  erratum. 
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SUMMARY 


The  analysis  of  the  stability  derivatives  of  low-aspect— ratio 
triangular  wings  at  subsonic  and  supersonic  speeds,  given  in  NACA  TN 
No.  Ih23,  is  extended  to  apply  to  triangular  wings  having  large  vertex 
angles  and  traveling  at  supersonic  speeds.  The  lift,  rolling  moment 
due  to  sideslip,  and  damping  in  roll  and  pitch  for  this  more  general 
case  have  been  treated  elsewhere  on  the  basis  of  the  theory  of  small 
disturbances.  The  surface  potentials  for  angle  of  attack  and  rolling 
taken  therefrom  are  used  to  obtain  the  lift  due  to  downward  acceleration, 
the  several  side— force  and  yawing-moment  derivatives  that  depend  on 
leading-edge  suction,  and  a  tentative  value  for  the  rolling  moment  due 
to  yawing.  All  the  known  stability  derivatives  of  the  triangular  wing 
at  supersonic  speeds,  regardless  of  source,  are  summarized  for  convenience 
and  presented  with  respect  to  both  body  axes  and  stability  axes.  The 
results  are  limited  to  Mach  numbers  for  which  the  triangular  wing  is 
contained  within  the  Mach  cone  from  its  vertex.  The  spanwise  variation 
of  Mach  number  in  the  case  of  yawing  is  neglected,  although  the  effect 
must  be  of  importance. 


INTRODUCTION 


An  earlier  investigation  (reference  l)  has  provided  theoretical 
stability  derivatives  of  low-aspect— ratio  wings  of  triangular  plan  form 
at  subsonic  and  supersonic  speeds.  The  restriction  to  lew  aspect  ratio 
was  a  consequence  of  the  limitations  of  the  theory.  Several  Investigators 
have  since  obtained  pressure  distributions  for  angle  of  attack,  rolling, 
pitching,  and  sideslip  at  supersonic  speeds  (references  2  to  6  and 
unpublished  analyses) ,  without  restriction  to  low  aspect  ratio.  These 
derivations  have  employed  variants  of  the  linear  theory  of  supersonic 
flew  end  have.  In  fact,  constituted  import-ant  steps  In  the  development 
of  the  theory. 

If  the  rotations  are  taken  about  the  vertex,  the  pressure  distribu¬ 
tion  for  each  motion  in  the  more  general  case  is  found  to  have  the  same 
shape  as  the  corresponding  low-aspect— ratio  approximation,  so  long  as 
the  triangular  wing  is  contained  within  the  Mach  cone  from  the  vertex. 
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The  magnitudes  differ  by  factors  which  are  functions  solely  of  the  ratio 
of  the  tangent  of  the  semivertex  angle  of  the  triangle  to  the  tangent  of 
the  Mach  angle.  The  same  similarity  exists  between  the  distributions  of 
surface  potential.  It  is  thus  relatively  simple  to  extend  most  of  the 
derivations  of  reference  1  to  remove  the  restriction  to  low  aspect  ratio 
for  supersonic  speeds.  Such  an  extension  is  made  in  the  present  paper. 

The  lift-curve  slope,  the  damping  in  roll  and  pitch,  and  (in  effect) 
the  rolling  moment  due  to  sideslip  have  been  evaluated  in  references  2 
to  6,  so  that  the  principal  contributions  of  the  present  paper  are  the 
several  side— force  and  yawing-moment  derivatives  and  a  tentative  value 
of  the  rolling  moment  due  to  yawing.  All  the  known  stability  derivatives 
of  the  triangular  wing  at  supersonic  speeds,  regardless  of  source,  are 
collected  herein  for  convenience  and  presented  with  respect  to  both  body 
axes  and  stability  axes.  Wings  with  dihedral  are  not  treated  (although 
they  were  included  in  reference  l)  and  the  results  are  limited  to  Mach 
numbers  for  which  the  wing  is  contained  within  the  Mach  cone  from  its 
vertex. 


SYMBOLS 


x,y,z 

u,v,w 

V 

M 

M' 

B 

a 

3 

e 

u 


rectangular  coordinates  (fig.  1) 

incremental  flight  velocities  along  x— ,  y— ,  and  z-axes, 
respectively,  figure  2;  induced  flow  velocities  along 
x—,  y— ,  and  z-axes  of  figure  1,  respectively 

angular  velocities  about  x— ,  y— ,  and  z-axes,  respectively, 
figure  2 

flight  speed 

stream  Mach  number  (V/Speed  of  sound) 


component  Mach  number  normal  to  wing  leading  edge 

cotangent  of  Mach  angle  (v/m5  -  l) 
angle  of  attack  (Flight  w/Y) 
angle  of  sideslip  (Flight  v/v) 
semivertex  angle  of  triangle 
Mach  angle  ^cot- \/m^  —  1^ 


MC 


\flTc‘ 
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AP 

P 

a 

b 


local  pressure  difference  between  lower  and  upper  surfaces  of 
airfoil,  positive  in  sense  of  a  lift 

density  of  air 

semividth  of  triangle  at  distance  x  from  vertex 
span  (base  of  triangle) 

root  chord  (height  of  triangle) 

(  2  /W2  o 

mean  aerodynamic  chord  I  3  s  -  /  (Local  chord)  dy  = 

\  s  U  0 

edge  slope  (—  =  —  = 

Vx  dx  k  2c  J 


aspect  ratio 


(?) 


S 

0 

Tj  =  COS' 

*  .  ,/T 

E*  (BC) 

F*  (BC) 

E"(BC)  = 

G(BC)  = 


area  of  triangle 
velocity  potential 


(H 


-i  z 
a 


-  B2C2 


complete  elliptic  integral  of  the  second  kind  vith 
modulus  k: 


2  \Jl  —  k^sin2z  dz^ 


complete  elliptic  integral  of  the  first  kind  vith 

dz 


modulus  k 


'(Si 


\  0  \Jl  —  k2sin2z 


E*(BC) 


1  -  B2C2 


(1  -  2B2C2)E'(BC)  +  B2C2f»(BC) 


H(BC)  =  3G(BC)  -  2E"(BC) 


cvj|  m 
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I(BC) 


2(1  -  B2C2) 


(2  -B2C2)E*(BC)  -  B2C2F,(BC) 

j(BC)  =  E"(BC)l(BC)\Jl  -  B2C2 

constant  defined  in  equation  (l6) 
yawing  moment 
lateral  force 


m 


"0 


N 

8 


suction  force  per  unit  length  of  edge 
.  Lift' 

lift  coefficient 


pitching-moment  coefficient 


rolling-moment  coefficient 


yawing-moment  coefficient 


lateral— force  coefficient 


profile-drag  coefficient 


induced  surface  velocity  normal  to  wing  leading  edge 
perpendicular  distance  of  point  (x,y)  from  wing  leading  edge 


eg 

Subscripts 


distance  of  center  of  gravity  forward  of 


B 

L 


right  edge 
left  edge 
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Whenever  a,  a,  q,  p,  |3,  and  r  are  used  as  subscripts,  a 
nondimens ional  derivative  is  indicated  and  this  derivative  is  the  slope 
through  zero.  For  example. 


A  dot  above  a  symbol  denotes  differentiation  with  respect  to  time. 
All  angles  are  measured  in  radians. 


ANALYSIS 

SCOPE 


The  stability  derivatives  of  triangular  wings  at  supersonic  speed 
that  have  been  treated  theoretically  herein  or  elsewhere  are  listed  in 
table  I,  together  with  the  expressions  that  have  been  found  for  them. 

All  the  derivations  make  use  of  body  axes.  The  derivations  that  follow 
give  the  values  with  reference  to  the  principal  body  axes  of  figure  2 
with  origin  at  the  aerodynamic  center  0,  o) .  Conversion  has  been 

made  to  the  system  of  stability  axes  shown  in  figure  3  with  origin  a 

distance  x„_  ahead  of  the  point.  Table  I  comprises  parallel  columns 
eg  3 

which  present  formulas  relative  to  both  systems.  The  expressions  are 
limited  to  Mach  numbers  for  which  the  triangle  is  contained  within  the 
Mach  cone  from  its  vertex. 


DERIVATIVES  CT  ,  CT . ,  and  C 

La'  La  ma 

The  pressure  distribution  on  a  thin  delta  wing  at  an  angle  of  attack 
in  a  supersonic  stream  has  been  obtained  in  references  2  to  4  by  the 
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linearized,  theory  without  restriction  on  the  vertex  angle  of  the  triangle. 
The  approximation  originally  given  for  the  slender  (low-aspect— ratio) 
triangle  (reference  7)  and  used  as  the  basis  for  reference  1  is  found  to 
apply  to  the  general  case  upon  division  by  a  constant  (an  elliptic  integral) 
that  depends  on  the  ratio  of  the  3emivertex  angle  to  the  Mach  angle.  That 
is. 


AP  _  _ taCa  ^ 

i  pV2  E»(BC)  ^a2  -  y2 

where  E*(BC)  is  the  complete  elliptic  integral  of  the  second  kind  with 
modulus 


k  =  \J  1  -  B^2 
=  \jl  -  (m2  -  ])  C2 


Thus  the  lift-curve  slope  for  the  more  general  case  is  the  value  given  by 
references  7  and  1  divided  by  E*(BC): 


rtA 

2E*(BC) 


=  |  AE"  (BC) 


(2) 


The  surface  potential  given  in  equation  (3)  of  reference  1  is  likewise 
extended  to  include  nonslender  triangles  at  supersonic  speeds  upon  division 
by  E*(BC).  The  revised  potential  is 


+  Vaa  sin  q 
E'(BC) 


E»(BC) 


(3) 


The  elliptic  integral  E*(BC)  depends  only  on  the  parameter  BC  *= 

(ratio  of  the  tangent  of  the  semi vertex  angle  of  the  delta  wing  to  the 


tan  e 
tan  u 
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tangent  of  the  Mach  angle)  and  la  therefore  a  constant  for  a  given  wing 


at  a  given  speed.  The  derivations  in  reference  1  for  CL> 

which  are  based  on  the  potential  0,  thus  merely  acquire  a 

E"(bc)  -  1 

E*  (BC) 

CL.  =  ~  Em(BC) 
a  2 

and  Cjj^, 

factor 

(b) 

c  ^  e"(bc) 

ma  1^ 

(5) 

DERIVATIVES  Cm  .  CT  ,  and  C, 

q  zp 

The  derivatives  C_,  ,  Ct  ,  and  C7  are  derived  in  reference  5* 

h  q 

With  respect  to  the  axes  of  figure  2 

°n  --^G(BC) 

q  16 

(6) 

CL  =  ^  H(BC) 

<1  2 

(7) 

c  =  -  I(BC) 

lV  32 

(8) 

where 

G(BC)  -  1  “  b2°2 

(1  -  SB^^E'fBC)  +  B^C%*(BC) 

(9) 

H(BC)  =  3G(BC)  -  2E"(BC) 

(10) 

.  .  2(1  -  B2C2) 

i(bc  )  - - - - - 

(2  -B2C2)E*(BC)  -BSc^^BC) 

(11) 
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and.  F*(BC)  and  E*(BC)  are  the  complete  elliptic  integrals  of  the  first 
and  second  kinds,  respectively,  with  modulus  k  =  ^1  —  B2c^. 


DERIVATIVE  C 

The  pressure  distribution  over  a  thin  delta  wing  in  yaw  (sideslip)  at 
an  angle  of  attack  at  supersonic  speed  has  been  obtained  in  reference  6  , 

and  unpublished  work.  If  the  angle  of  yaw  is  assumed  to  be  small  fp  «  l' 

the  rolling-moment  coefficient  can  be  expressed  in  the  approximate  form  ' 


Cz  ~_2|&e"(BC) 

Thus,  the  derivative  with  respect  to  3  is 


C,  =  —  E" (BC)  (12) 

3 

An  alternative  derivation  based  on  the  surface  potential,  equation  (3), 
for  the  unyawed  wing  will  be  given  because  the  method  provides  the  starting 

point  for  a  derivation  of  C7  ,  C_  ,  and  Cv  . 

4r  ”r  Yr 

The  potential  for  the  disturbance  velocity  may  be  expressed  relative 
to  axes  alined  with  the  stream  (wind  axes)  or  with  respect  to  axes  that 
yaw  with  the  body  (body  axes).  For  small  angles  of  yaw  (p  «  the 

linearized  equation  for  the  potential  has  the  same  form  relative  to  either 
system  of  axes.  The  potential  is  determined  by  the  normal  velocity  of 
points  of  the  surface  and  by  the  orientation  of  the  surface;  for 
negligible  thickness,  this  normal  velocity  is  Just  aV  for  all  angles 
of  yaw.  The  potential  expressed  relative  to  wind  axes  thus  varies  as 
the  wing  yaws  relative  to  these  axes.  The  potential  expressed  relative 
to  body  axes  is  constant  for  3mall  yaw  because  the  orientation  of  the  wing 
relative  to  the  axe 3  does  not  change. 

For  wind  axes,  Bernoulli*s  law  has  the  form 


AP  =  2pV  M 

ox 

and  the  change  in  the  pressure  distribution  with  yaw  results  from  the  chan^ 
in  the  potential  function  wibh  yaw.  For  body  axes  with  small  yaw, 
Bernoulli* s  law  has  the  approximate  form 
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AP  =  2pV(  M  -  3  ^ 
\ox  oy> 


(13) 


and  the  change  in  pressure  distribution  with  yaw  results  from  the 
term 


-3  ^  since  0  does  not  change. 


In  reference  1  in  the  section  entitled  "DERIVATIVE  C7  ,"  the 

3 

derivation  employs  body  sixes  and  equation  (13)  of  the  present  paper. 

The  surface  potential  used  (equation  (3)  of  reference  1)  is  the 
approximation  for  narrow  vertex  angle.  Equation  (3)  herein  for  a 
general  vertex  angle  may  be  used  instead.  Equation  (3)  herein  differs 
only  in  the  factor  l/E*(BC),  whence  the  earlier  expression  for  C7 

P 

(equation  (19),  reference  1,  with  r  =  0°)  acquires  this  factor  to  agree 
with  equation  (12). 


DERIVATIVE  C 

*r 

The  foregoing  discussion  of  the  triangular  wing  in  yaw  (sideslip) 
may  be  extended  to  provide  a  preliminary  treatment  of  the  case  of  a 
small  angular  velocity  of  yaw  r.  The  corresponding  extension  for 
narrow  vertex  angle  is  made  in  reference  1.  The  treatment  is  general¬ 
ized  to  an  arbitrary  vertex  angle  for  supersonic  speeds, as  before,  by 
using  equation  (3)  herein  for  the  surface  potential.  Two  changes  then 
appear  in  the  pressure  equation,  equation  (20),  of  reference  1.  The 
right-hand  side  is  divided  by  E*(BC),  and  the  term  aC  =  xC2  must  be 
retained,  since  C“  is  no  longer  small  compared  with  unity  (C  =  tangent 
of  semivertex  angle).  With  these  changes,  the  derivation  leads  to 

C,  =  nix  (—  +  ~ |E"(BC)  (14) 

lr  \9A  16  J 

In  the  derivation  of  equation  (14),  the  spanwise  variation  in  local 
Mach  number  caused  by  yawing  is  not  taken  into  account  although  the 
variation  in  forward  speed  is  taken  into  account.  The  surface  potential 
that  is  used,  equation  (3),  satisfies  the  linearized  equation  for  a  flow 
of  uniform  Mach  number.  This  potential  is  inadequate  to  describe  the 
compressibility  effects  associated  with  a  spanwise  variation  of  Mach 
number. 

Thus,  consider  a  high— aspect— ratio  rectangular  wing  with  tips  cut 
off  along  the  Mach  lines.  In  straight  flight  the  Ackeret  theory  can  be 
applied.  The  pressure  difference  is  given  by 
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AP  = 


2oM2(Speed  of  sound)2 

a - 


(15) 


In  yawing  flight  the  forward  velocity  varies  linearly  along  the  span.  If 
the  rate  of  yaw  is  made  sufficiently  low,  the  variation  from  wing  tip  to 
wing  tip  can  be  made  so  small  that  the  flew  is  still  nearly  two-dimensional, 
at  any  point.  Thus  the  Ackeret  theory  is  still  applicable  if  the  local 
Mach  number  is  used  at  each  spanwise  station. 

The  variation  in  pressure  with  local  Mach  number  can  be  obtained  from 
equation  (15).  As  the  Mach  number  is  increased,  the  pressure  decreases 
from  infinity  at  M  =  1  to  a  minimum  at  M  =  1.4  and  then  increases  again. 
Thus  below  Mach  number  1.4  the  faster  moving  sections  of  the  yawing  wing 
have  the  lesser  lift.  This  result  is  contrary  to  subsonic  behavior  and  to 
that  which  would  be  predicted  if  the  spanwise  variation  of  Mach  number  were 
neglected.  Thus  the  spanwise  variation  of  the  compressibility  effect  causes 
a  reversal  of  the  sign  of  the  rolling  moment  due  to  yawing  for  rectangular 
wings  at  Mach  numbers  between  1  and  1.4,  and  at  M  =  1.4  the  moment  is 
zero.  (This  result  refers  to  yawing  in  a  system  of  stability  axes,  fig.  3. 
For  body  axes,  fig.  2,  the  effect  is  similar  but  the  reversal  extends 
to  M  =  °°.  ) 

A  yawing  triangular  wing  may  be  expected  likewise  to  Bhow  an  effect 
of  the  spanwise  variation  in  Mach  number.  If  the  triangle  is  contained 
within  the  Mach  cone  from  its  vertex  (the  only  case  considered  in  this 
paper),  however,  the  effect  should  be  very  much  less  than  for  the  rec¬ 
tangular  wing.  In  particular,  where  the  predicted  effect  for  the  rec¬ 
tangular  wing  is  a  reversal  of  the  sign  of  the  rolling  moment,  the  effect 
for  the  triangular  wing  is  expected  to  be  merely  a  change  in  the  magnitude. 

A  reversal  in  sign  is  not  expected  until  the  edges  of  the  triangle  protrude 
from  the  Mach  cone.  This  behavior  Is  inferred  from  the  fact  that  the 
analyses  of  references  2  to  7  show  many  subsonic  characteristics  for 
triangles  within  the  Mach  cone  and  a  marked  change  in  characteristics  for 
triangles  with  side  edges  outside  the  Mach  cone. 


DERIVATIVES  Cv  and  Cn 

h  H? 

Extensive  changes  are  necessary  to  generalize  the  treatment  of  CY 

P 

and  in  reference  1  to  arbitrary  vertex  angles  for  supersonic  speeds; 

therefore,  the  revised  derivation  Is  given  in  detail. 

The  derivatives  Cy  and  Cn  relative  to  body  axes  for  a  very 

p  P 

thin  delta  wing  without  dihedral  arise  entirely  from  suction  on  the  wing 
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side  edges.  Consider  a  condition  for  which,  the  induced  velocity  normal 
to  the  edge  is  of  the  form 


(16) 


in  the  immediate  neighborhobd  of  the  edge,  where  s  is  the  perpendicular 
distance  from  the  edge  and  K  is  a  constant.  Reference  3  points  out 
that  for  such  a  flow  there  is  a  suction  force  per  unit  length  of  edge. 


f  = 


(17) 


so  long  as  the  delta  wing  does  not  protrude  from  the  Mach  cone  from  its 
vertex.  In  equation  (17),  M*  is  the  Mach  number  of  the  component  of 
the  stream  flow  normal  to  the  leading  edge.  The  radical  \J  1  —  M'2  is 
the  Prandt 1-Glaue rt  compressibility  factor  for  the  normal  component  of 
flow.  Equation  (17)  is  limited  to  real  values  of  the  radical  by  the 
condition  expressed  for  the  Mach  cone. 

For  the  delta  wing  in  rolling  motion  the  induced  velocity  component 
has  been  obtained  in  reference  5  as 


Angle  of  attack  gives  the  additional  contribution  (reference  2) 

dVC2 

u  =  ±  —  —  — — 

y  E*(BC)^C2  r (zy 


The  total  induced  velocity  on  the  upper  surface  is  thus  the  sum  of 
and  u^  with  the  plus  sign 


Very  near  the  side  edge  this  velocity  is  approximately 


12 
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u 


«) 


-22—  i  IS  i(bc) 

E'(BC)  2 


where  the  plus  sign  refers  to  the  right  edge  and  the  minus  sign  to  the  left 
edge. 


If  a  similar  calculation  is  made  for  v 


it  is  found  that  as 


the  side  edge  is  approached  the  resultant  Induced  velocity  +  v2 

becomes  normal  to  the  edge.  Thus  the  normal  velocity  near  the  edge  is 


v 


H  " 


\Jl  +  cg 

c 


The  perpendicular  distance  of  point  (x,y)  from  the  side  edge  is 


s 


<c  -  HD 

\Ji  + 


The  resultant  induced  velocity  very  near  the  edge  may  therefore  be 
expressed  approximately  as 


aV  .  l(BC)pCx 

..  _  T  — 

E'(BC)  "  2 


which  is  of  the  form  of  equation  (16).  The  suction  force  per  unit  length 
of  edge  is  from  equation  (17)  thus 


f  = 


n 

2 


l(BC)]gp2C2x2 

4 


+  !(BC)aVpCx 
“  E*(BC) 


\j(l  +  C2)(l  -  M*2) 


where  the  plus  sign  refers  to  the  right  edge  and  the  minus  sign  refers  to 
the  left  edge .  The  factor  \j ( 1  +  C2)(l  —  M*2)  can  be  reduced  to 
sjl  —  B^C  ,  where  B2  =  M2  -  1. 

The  lateral  component  of  this  suction  force  is  given  by 
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’Jo  (fE-fL) 


K  po  l(BC)  \J1  —  b2c£ 

=  I  pC2c3aVp  — - — E - 

3  E'(BC) 


1  2 

The  lateral— force  coefficient  is  formed,  by  division  by  —  pV  S,  and 

2 

the  derivative  with  respect  to  pb/2V  is  the  stability  derivative  Cy  . 
It  is  p 


2na  I(BC)  -  B2C£ 


'Yn  "  3 


E*(BC) 


The  yawing  moment  of  the  leading-edge  suction  about  the  vertex  of 
the  triangle  is 


bo  ■  -L  C1 


fE-fL)  xJTT?  dxJl  +  C2 


=  -  f  pC2c^aVp(l  +  C2) 


2  i(bc)  yi  —  b2cc 
e*(bc) 


The  moment  about  the  reference  point 


(fc,0,o) 


N  =  +  —  cY 

0  3 


*  p  h  /,  I(BC)  Jl  -  B2C£ 

= - pC^c  <xVp(l  +  9C2}  — - — - - - 

36  E*(BC) 


1  2 

The  yawing-moment  coefficient  is  formed  by  division  by  —  pVSb,  and 

the  derivative  with  respect  to  pb/2V  is  the  stability  derivative  Cr 
It  is 


<v ,  _  Jl.  *  ms)A  ~.i!f 

“P  \9A  16/  E'(BC) 
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DERIVATIVES  CY  ,  C  , 

V  n3 


C„  ,  and  C 
*r 


According  to  the  discussion  on  C,  ,  a  small  angle  of  yaw  or 

/  \  *3 

sideslip  (3  «  j-jj  does  not  alter  the  surface  potential  expressed 

relative  to  body  axes.  As  a  consequence  the  induced  velocity  distribu¬ 
tion  is  unchanged.  Thus  the  initially  symmetric  distribution  of  leading- 
edge  suction  persists  in  sideslip,  and  the  derivatives  Cy  and  Cn 
are  zero.  3  3 


The  surface  potential  relative  to  body  axes  is  likewise  unaltered  in 
the  first  approximation  by  a  small  angular  velocity  of  yaw.  Accordingly, 
insofar  as  the  pressure  forces  are  concerned,  the  derivatives  CY  and  C 

*r  ¥ 

are  zero.  Subsonic  experience,  however,  suggests  an  appreciable 

Cn  -derivative  (damping  in  yaw)  from  profile  drag.  This  damping  derivative 
¥ 

has  been  evaluated  in  reference  1  as 


<V--Cd/^  — 

^  V  9Ae 


RESULTS  AND  DISCUSSION 

The  formulas  that  have  been  obtained  for  the  various  stability 
derivatives  are  collected  In  table  I.  Derivatives  obtained  elsewhere 
are  Included  for  completeness,  and  the  source  is  Indicated  In  each 
instance.  Expressions  are  given  for  two  systems  of  coordinate  axes. 

In  the  first  column  are  shown  the  derivatives  relative  to  the  principal 


body  axes  of  figure  2  with  origin  a  distance 


from  the  vertex  of  the 


triangle.  In  the  second  column  are  shown  the  results  relative  to 

stability  axes  with  origin  a  distance  x  ahead  of  the  ~c  poll 

eg  3 


point . 


The  relationship  between  the  two  systems  of  axes  is  shown  in  figure  3. 
Equations  for  transforming  from  body  axes  to  stability  axes  are  given 
in  reference  8;  the  shift  In  origin  results  in  additional  terms. 


In  the  transformation  of  the  present  results  from  principal  body 
axes  to  stability  axes  terms  of  order  A^/l6  and  the  more  important 
terms  of  order  a?  are  retained  (see  footnote,  table  I),  whereas  in 
reference  1  such  terms  are  dropped  as  a.  consequence  of  the  narrow  vertex- 
angle  approximation. 
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These  results  for  an  arbitrary  vertex  angle  may  be  compared  with  the 
asymptotic  values  for  the  case  of  vertex  angle  approaching  zero  given  in 
reference  1.  The  present  results  for  principal  axes  are  found  to  differ 
from  the  asymptotic  values  (except  for  small  terms  in  P?  and  aP)  only 
In  the  acquisition  of  certain  factors  which  are  functions  of  BC.  Thus 

the  asymptotic  values  for  CT  >  Ct  . ,  C™  ,  C7  ,  and  C7  are  multiplied 

■4x  q  4  3  ir 


by  E"(BC);  Cm  is  multiplied  by  G(BC);  Ch  is  multiplied  by  H(BC); 

Q  q 

C7  is  multiplied  by  l(BC);  and  C_  and  Cy  are  multiplied 


I  (BC)  \ll  -  B2C2 


by  ■  - - - -  =  J(BC).  The  parameter  BC  =  — == — *  is  the  ratio  of 

E'(BC)  tan  ^ 

tangent  of  the  semivertex  angle  of  the  triangle  to  the  tangent  of  the 
Mach  angle.  BC  approaches  zero,  therefore,  as  the  vertex  angle  approaches 
zero.  The  several  functions  E" (BC),  .  .  .,  J(BC)  all  approach  unity 
as  BC  approaches  zero,  and  thus  the  derivatives  obtained  herein  approach 
the  asymptotic  values  of  reference  1  as  the  vertex  angle  goes  to  zero. 


The  variation  of  these  stability  derivatives  with  Mach  number  is 
contained  entirely  in  the  factors  E"(BC),  ...»  J(BC).  The  five 

factors  are  plotted  against  BC  =  ■L-^n--£,  the  ratio  of  the  tangent  of 

tan  p 

the  semi vertex  angle  to  the  tangent  of  the  Mach  angle,  in  figure  4. 


The  derivatives  apply  to  a  wing  of  triangular  plan  form  and  zero 
thickness.  The  calculations  are  based  on  the  assumption  of  potential 
flow  with  small  disturbances,  except  in  the  case  of  the  derivative  Cn  , 

in  which  shin  friction  is  considered.  The  predicted  infinite  negative 
pressure  acting  on  an  edge  of  zero  thickness  to  yield  a  finite  suction 
force  is,  of  course,  a  mathematical  idealization.  (The  local  violation 
of  the  assumption  of  small  disturbances  is  not  serious.)  Subsonic 
experience  indicates  that  with  a  suitably  rounded  edge  a  considerable 
leading— edge  suction  force  may  be  realized  in  practice,  with  the 
theoretical  value  an  upper  limit.  On  the  other  hand,  a  sharp  leading  edge 
Is  known  to  cause  loss  of  the  leading-edge  suction.  The  requirements  of 
extreme  thinness  and  a  rounded  leading  edge  (that  is,  appreciable  radius 
of  curvature)  are  evidently  in  conflict.  Thus  the  degree  of  applicability 
of  the  yawing-moment  and  lateral -force  derivatives  to  actual  triangular 
wings  is  uncertain.  A  further  limitation  on  validity,  already  elaborated 

on  in  the  section  on  C?  ,  exists  also  for  the  derivatives  with  respect 

Lr 

to  yawing  velocity.  The  analysis  neglects  the  spanwise  variation  in 
Mach  number  caused  by  the  yawing  (but  not  the  spanwise  variation  in 
velocity).  The  result  is  an  error  in  the  magnitude  of  the  yawing 


1 6 
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derivatives  that  is  expected  to  vary  from  zero  for  BC— >0  to  an  important 
amount  for  BC — >1. 

The  potential  0  satisfies  the  linearized  equation  of  motion  for 
the  steady  state  but  not  the  more  general  linearized  equation  for  unsteady 
motion.  This  circumstance  implies  that  the  present  expressions  for  the 
stability  derivatives  are  suitable  only  for  steady  motions,  motions  with 
small  accelerations,  or  sinuous  motions  of  low  frequency.  This  limitation 
is  accepted  in  all  stability  work  and  may  become  serious  only  in  cases  of 
high-frequency  oscillations  such  as  flutter. 
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TABLE  I.-  STABILITY  DERIVATIVES  OF  THIN  TRIANGULAR  WINGS  AT  SUPERSONIC  SPEEDS 


Principal  body  axes 


Stability  axes  . 

origin  at  distance  xfig  ahead  of  |c  pointj 
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Figure  1.-  Axes  and  notation  used  in  analysis. 


Figure  2.-  Velocities,  forces,  and  moments  relative  to 

2 

principal  axes  with  origin  at  — c . 


Figure  5.-  Velocities,  forces,  and  moments  relative  to 

2 

stability  axes  with  origin  at  — c  -  xC2.  Principal 
axes  of  figure  2  dotted  in  for  comparison. 
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Mach  I  me 


E  (BO 
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H(BC) 


Figure  <4.  __  E/hpt/c  m  tegra /  factors  of  the 
stability  derivatives  that  determ  me 
their  variation  with  Mach  number. 

(  See  fable  X. ) 


